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The temporal response of resonances in nanoplasmonic structures typically converts an incoming 
few-cycle field into a much longer near-field at the spot where non-linear physical phenomena in- 
cluding electron emission, recollision and high-harmonic generation can take place. We show that 
for practically useful structures pulse shaping of the incoming pulse can be used to synthesize the 
plasmon-enhanced field and enable single-cycle driven nonlinear physical phenomena. Our method 
is demonstrated for the generation of an isolated attosecond pulse by plasmon-enhanced high har- 
monic generation. We furthermore show that optimal control techniques can be used even if the 
response of the plasmonic structure is not known a priori. 



PACS numbers: 42.65.Ky, 78.67.Bf, 73.20.Mf 

In the last decade, there has been rapid progress 
in the production of ultrashort light pulses with con- 
trolled waveforms and durations down to below 100 at- 
toseconds [1-3]. These light sources enable the con- 
trol and tracing of electron dynamics in atoms [4- 
6], molecules [7, 8], and solids [9, 10] on their natural 
timescale. The realization of a similar level of control 
of the electron motion in nanocircuits has the poten- 
tial to revolutionize modern electronics [1]. Lightwave 
controlled nanocircuits (lightwave nanoelectronics) may 
reach petahertz operation frequencies and might remove 
the bottleneck in conventional communication technol- 
ogy by enabling all-optical information processing and 
communication. The key to lightwave nanoelectronics 
is the control of electron dynamics in nanostructured 
materials on sub-cycle timescales. Progress has very re- 
cently been made in the control of electron dynamics in 
nanoparticles [11], nanotips [12] and nanojunctions [13] 
with carrier-envelope phase (CEP) stabilized few-cycle 
pulses. The control of few-cycle waveforms by the CEP, 
however, gives only a very limited degree of control over 
the electron dynamics, which can be significantly im- 
proved with ultrabroadband lightwave synthesis permit- 
ting to sculpt the electric field of a laser pulse with at- 
tosecond precision [14]. Despite its importance for light- 
wave nanoelectronics, this approach has not yet been im- 
plemented for the shaping of plasmonic near- fields. 

The application of few-cycle pulses to a resonant 
nanostructure typically results not only in the desired 
field-enhancement effect which may be utilized for non- 
linear applications but also in a temporally longer near- 
field evolution [15]. It is desirable, however, to reach 
single-cycle near-field profiles enabling processes such as 
a well controlled asymmetry in the direction of electronic 
currents, the confinement of electron emission and ac- 
celeration to a single cycle, and the generation of iso- 
lated attosecond extreme ultraviolet (XUV) pulses. Fur- 



thermore, if CEP-control of strong-field processes that 
are enabled by the near-fields is required, the amplitude 
of those CEP-effects typically scale inverse exponentially 
with the near-field duration and their relevance is thus 
limited to the few-cycle regime [16]. 

In this Letter, we will present an approach employing 
pulse shaping techniques to form plasmonic near-fields, 
which enable nonlinear phenomena that are induced by 
just a single cycle of the near-field. While this approach 
is general and not limited to a certain nonlinear pro- 
cess, we will discuss its implementation for the generation 
of isolated attosecond pulses via nanoplasmonic field en- 
hancement. Plasmon- assisted high-harmonic generation 
(HHG) was studied experimentally [17, 18] and theoret- 
ically [15, 19]. Here, the strong field enhancement ob- 
tained from plasmon resonances in metallic nanostruc- 
tures is used to locally enhance the electric field strength 
to the levels required for HHG. As the local intensity 
can be enhanced by more than 4 orders of magnitude, 
this drastically lowers the required intensity of the driv- 
ing laser pulses and enables HHG for repetition rates in 
the MHz range. The approaches, however, do not cur- 
rently produce isolated attosecond pulses, which are a 
critical ingredient for attosecond spectroscopy, since the 
resulting near-field oscillations are too long. It should 
be noted that the experimental studies in [17] were re- 
cently challenged by Sivis et al, who only observed XUV 
fluorescence from bow-tie nanoantennas for similar con- 
ditions [20] . Kim et at. show in recent work that the local 
rare gas pressures can be increased sufficiently to reach 
the regime of coherent XUV generation by HHG [21]. 
Further studies will be required to clarify the importance 
of coherent versus incoherent XUV emission from such 
nanostructures. Since HHG is only used as an example 
application for nanoplasmonic near-field synthesis in our 
studies, we will neglect incoherent processes in the dis- 
cussion below. 
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In order to support the generation of isolated attosec- 
ond pulses, the high-harmonic generation process has 
to be gated so that it only occurs within a short win- 
dow of time. Husakou et at. have recently shown that 
one common approach, polarization gating [22-24], can 
be transferred to the nanoscale for plasmonic structures 
that replicate the polarization properties of the incom- 
ing pulse [19]. We here propose a more general approach 
where the incoming pulse is shaped such that the locally 
enhanced field at the nanostructure achieves temporal 
gating of the HHG process. This approach is quite inde- 
pendent of the specific properties of the plasmonic struc- 
ture. In particular, we find that even when the lifetime 
of the plasmon resonance is much larger than the cycle 
time of the IR field, the local response can be shaped 
in such a way as to allow for the generation of isolated 
attosecond pulses with large contrast. 

As an example, we study one of the structures inves- 
tigated in [15]: two gold ellipsoids with major and mi- 
nor axes of 100 nm and 16.7 nm, respectively. The major 
axes of both ellipsoids are aligned along the same axis, 
with a gap of 5nm between them. This geometry cre- 
ates a local "hot spot" between the ellipsoids, where the 
field enhancement is expected to be maximal. We focus 
on the technique of amplitude gating [25, 26], where the 
generating IR field only becomes strong enough to gen- 
erate high-energy harmonics in the cutoff region during 
a short time window. An isolated attosecond pulse can 
then be generated by spectrally filtering the generated 
high harmonic radiation. 

The gold ellipsoids are described in the limit of linear 
dielectric response, with a frequency-dependent dielectric 
function e r (uj) taken from the experimental (bulk) values 
in [27]. The strongly nonlinear response of the atomic 
gas used as the HHG medium is small enough that it is 
neglected when solving the Maxwell equations. In ad- 
dition, since the generation of high harmonic photons is 
confined to scales far below the driving wavelength, co- 
propagation and phase matching of the harmonics and 
driving pulse do not have to be explicitly treated. Be- 
cause of the linearity of the Maxwell equations in this 
approximation, the plasmon response to an ultrashort 
(broad-band) incoming pulse can be calculated as the 
superposition of fixed-frequency components with appro- 
priate amplitudes. The spatiotemporal electric field dis- 
tribution is then given by 

E(f,t) = ^c nJ B„(r)e- iw "* +c.c, (1) 

n 

where E n (f) is the (complex- valued) spatial response for 
incoming mode n at frequency cj n , and c n are the complex 
amplitudes determining the temporal shape of the incom- 
ing pulse. We choose all incoming modes to be plane 
waves polarized along the axis connecting the ellipsoids 
(the z-axis) and propagating along the same orthogonal 




FIG. 1. Frequency-dependent intensity enhancement (blue) 
and phase shift (red) at the center of the structure. The insets 
show the spatial dependence of enhancement and phase shift 
for three frequencies: In the quasistatic regime, close to the 
resonance peak, and on the blue side of the resonance. 

direction (x-axis). Thus, E n (r) « ze lWnX ^ c for large neg- 
ative x. Using plane waves is equivalent to assuming that 
the focus spot size is large compared to the extension of 
the system. The response E n (r) for each fixed-frequency 
component is obtained by solving a numerical scatter- 
ing problem using the SCUFF-EM package [28, 29], a free, 
open-source software implementation of the boundary- 
element method (BEM) of classical electromagnetic scat- 
tering [30]. The BEM exploits known Maxwell solutions 
to express the fields inside and outside homogeneous ma- 
terial bodies in terms of effective surface currents flow- 
ing on the body surfaces. This has the advantage that 
we need only discretize surfaces, not volumes, yielding a 
computationally efficient approach. 

The frequency components used to construct (synthe- 
size) the pulses are integer multiples of uq = 0.031 eV 
(giving a frequency comb). The incoming field is thus 
a periodic pulse train with period T = 2tt/ujo « 133 fs, 
which is long enough to ensure that the response of the 
system to a single pulse has decayed by the time the next 
one arrives. We checked that using half the period leaves 
the results essentially unchanged. The periodicity thus 
does not influence our conclusions. 

We start by discussing the response E n (f), which in 
principle contains all information about the system (in 
the limit of linear response). Fig. 1 shows that there is a 
dominant plasmon resonance peak at a photon energy of 
around 1.443 eV (A w 859 nm), with a lifetime of ~7.5fs. 
The maximum intensity enhancement at the center be- 
tween the ellipsoids is « 3 • 10 4 . As the resonance is 
crossed, the phase jumps by tt as expected (illustrated 
by the dominant z-component of the field). When driv- 
ing with a short (broad-band) pulse, the response will 
thus be in phase for the frequency components on the 
red side of the resonance, but out of phase on the blue 
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side. This implies that a Fourier- limited few-cycle in- 
coming pulse will produce a distorted and prolonged re- 
sponse in the center of the structure, essentially because 
the plasmon resonance is excited and keeps oscillating 
even after the driving pulse is over [see Fig. 2(a)]. It is 
this distortion that we wish to compensate in the follow- 
ing. The insets in Fig. 1 show the spatial dependence 
of the plasmon response of the structure. The field is 
strongly enhanced at the sharp tips of the ellipsoids, with 
the largest enhancement in the narrow gap at the center. 
HHG will therefore be dominantly produced in this "hot 
spot". Crucially, although not surprisingly, the phase of 
the response is spatially uniform within this hot spot for 
all relevant frequencies. Fourier synthesis of the incom- 
ing pulse will thus produce a uniform near-field temporal 
response, enabling control of strong-field processes. 

To produce isolated attosecond pulses with amplitude 
gating, the driving laser pulse must only reach a suffi- 
ciently large electric field strength for a single cycle, and 
with the correct carrier-envelope-phase [1, 25, 26, 31, 32]. 
Once the frequency-dependent response is known, it is 
in principle possible to choose the incoming pulse (de- 
scribed by the amplitudes c n in Eq. 1) such that an 
arbitrary pulse shape is created in the hot spot where 
HHG can then take place. This could be achieved with 
standard pulse shapers that allow control over amplitude 
and phase of the separate frequency components of an 
incoming broadband pulse [33, 34]. The "conventional" 
pulse shape desired in the hot spot would be a few-cycle 
Gaussian pulse as commonly used in amplitude-gated iso- 
lated attosecond pulse generation. The field response for 
an incoming Fourier-limited Gaussian pulse and for such 
an adjusted pulse is demonstrated in Fig. 2(a,b). How- 
ever, this method of field synthesis has some drawbacks: 
Since the plasmon resonance is typically more narrow- 
band than the bandwidth of a few-cycle pulse (i.e., the 
lifetime of the resonance is longer than the desired pulse) , 
the incoming pulse must have a "hole" in the spectral 
distribution at the frequencies that are most strongly en- 
hanced, such that the enhanced pulse becomes Gaussian. 
This most straightforward approach will therefore limit 
the achieved field amplification. In addition, it requires 
that the structure is perfectly characterized, which in re- 
ality will not always be the case. A more systematic 
approach that circumvents these limitations is therefore 
preferable. For this we turn to the well-developed tools 
of optimal control to automatically generate the field in 
the center of the structure so that an isolated attosecond 
pulse is produced with the desired parameters. 

We assume the following setup: A laser system pro- 
duces an incoming broadband Gaussian pulse with a 
(transform-limited) peak intensity of 4-10 11 W/cm 2 , with 
a central energy of 1.448 eV (A « 856 nm), close to the 
plasmon resonance, and a FWHM bandwidth of 0.55 eV, 
corresponding to a FWHM duration of 3.33 fs. This pulse 
is sent through a pulse shaper which can separately con- 
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TABLE I. Properties of the attosecond pulses produced by 
HHG with the plasmon-enhanced field for different incoming 
pulses and filter settings. In all cases, the total intensity and 
CEP have been optimized. "Noise" is the relative power in side 
pulses, i.e., 1 — P as /Ptot- Filter setting "none" corresponds 
to using an unfiltered few-cycle incoming pulse, as in [15], 
while "flat phase" means compensating (only) the phase of 
the incoming pulse so the near field has a flat spectral phase. 
For these settings, the power in the incoming pulse does not 
influence the pulse shape. For filter setting "Gaussian", the 
amplitudes and phases are adjusted to produce a transform- 
limited Gaussian near field. The flanks of the peak are limited 
by the available spectral power and are thus not be perfectly 
Gaussian. The last two lines for each intensity show the fully 
optimized results and the results if the optimization is stopped 
after 10 4 function evaluations. 

trol the phases and magnitudes of the frequency compo- 
nents. The electric field in the structure is then 

E(r, t) = J2 c„/„£„(r>-^' + c.c. , (2) 

n 

where f n are complex numbers describing the setting of 
the pulse shaper (0 < \f n \ < 1). This field generates 
high harmonics in the hot spot of the two ellipsoids, 
which are then sent through a spectral filter centered 
at 100 eV with 10 eV FWHM bandwidth. This filtered 
XUV pulse is characterized by, e.g., attosecond streak- 
ing [35, 36], or alternatively a recently proposed method 
which only requires the HHG spectrum and thus seems 
particularly promising in light of its reduced measure- 
ment requirements [37]. To model this experimental 
setup, we simulate HHG using the strong-field approx- 
imation (SFA) [38, 39]. As we focus here on the con- 
trollability of the near field, we do not currently include 
additional effects induced by the presence of the metal 
surface [40]. 

The attosecond pulse generation is optimized through 
a feedback loop. The target function to be maximized 
determines the pulse properties that are optimized for. 
We here focus on two properties: i) The contrast, i.e., 
most of the power should be concentrated in a single 
attosecond pulse, and ii) the power in that pulse. The 
target function is F = P as /Ptot + r logP as , where P tot is 
the total power after filtering, P as is the power in the 
strongest attosecond pulse, and r is a coefficient cho- 
sen to produce the desired balance between pulse power 
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FIG. 2. Different synthesized pulses. In each row from left to right: Filter settings f n and electric field amplitudes E n (absolute 
value and phase are plotted separately), plasmon-enhanced electric field in the hot spot and filtered high harmonic radiation 
generated by this pulse (inset: incoming pulse), and time- frequency distribution of unflltered high harmonic radiation. The 
incoming pulse has a maximum intensity of 4 • 10 11 W/cm 2 before filtering. The rows, from top to bottom, show: (a) Fourier- 
limited incoming pulse with only amplitude and CEP adjusted to optimize single as pulse production in the hot spot, (b) 
Manually optimized pulse chosen to produce a few-cycle Gaussian pulse in the hot spot as well as possible, amplitude and CEP 
optimized, (c) Pulse obtained after numerical optimization (see text). The numerically optimized pulse features a phase jump 
before the "spike" producing the attosecond pulse, effectively suppressing HHG from the previous half-cycle of the field. 



and selectivity. We here choose r = 0.01, with powers 
in atomic units. This function is maximized by adjust- 
ing the coefficients f n in Eq. 2, using the Subplex al- 
gorithm as implemented in the NLopt package [41, 42]. 
The pulse and filter are described by 79 frequency com- 
ponents uj n = nujo with n G [8,86]. The optimization 
converges after around 10 5 function evaluations. How- 
ever, as the convergence is faster in the beginning and 
then slows down, the pulses obtained after 10 4 function 
evaluations are almost as good (cf. Table I) . In an actual 
experiment, the time needed to optimize the pulse may 
be under a few tens of seconds, if an all-optical character- 
ization of the pulses [43] is applied (limited by the rate at 
which HHG spectra can be acquired and analyzed, which 
can reach the kHz regime). 

We now turn to the results of the optimization proce- 
dure, shown in Fig. 2(c) and Table I, and general con- 
clusions that can be drawn. The exact pulse that is pro- 
duced varies depending on the parameters, but there are 



general features that we have found in all cases: The am- 
plitude gating (just one half-cycle with enough intensity 
to produce the desired harmonic frequencies) is achieved 
not by creating a few-cycle pulse, but through a relatively 
long pulse with a single superposed "spike". If we express 
the electric field as E(t) — f(t)exp(icj c t) + c.c, where 
uj c is the carrier or central frequency and f(t) is a slowly 
varying (complex) amplitude, the main spike can be char- 
acterized by a minimum and a phase jump of ~ tt in f(t) 
(similar to a zero crossing if f(t) was purely real), which 
creates a field minimum followed by a single maximum, 
cf. Fig. 2(c). The phase jump effectively suppresses HHG 
in the half-cycle before the pronounced spike that gen- 
erates the main attosecond pulse, which maximizes the 
single attosecond pulse generation efficiency, and espe- 
cially the contrast. This feature can conceptually be un- 
derstood as arising from the superposition of two pulses: 
One long background pulse, and a few-cycle pulse which 
is offset in central frequency and chirped. These two 
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pulses switch from destructive to constructive interfer- 
ence within the width of the few-cycle pulse. In this way, 
the few-cycle pulse can be much broader in frequency 
than the plasmon resonance and "boost" its intensity for 
HHG by constructive interference with a narrow-band 
(long) pulse. The long pulse by itself does not reach 
the intensity to produce high harmonic radiation at the 
100 eV desired for the attosecond pulse, so its longer du- 
ration does not influence the results. 

To summarize, we have introduced a general approach 
to synthesize few-cycle nanoplasmonic near-fields and 
demonstrated its potential for the generation of isolated 
attosecond pulses. We have shown that pulse shapes can 
be engineered to produce isolated attosecond pulses even 
in situations where the plasmon resonance is too long- 
lived to support generation of a strongly enhanced few- 
cycle near field, and that these pulse shapes can be man- 
ually chosen or automatically generated by an optimal 
control approach. 
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